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ABSTRACT 

The X-shooter instrument on the VLT was used to obtain spectra of seven moderate-redshift quasars simulta- 
neously covering the spectral range ^3000 A to 2.5 ^m. At z ~ 1 .5, most of the prominent broad emission lines 
in the ultraviolet to optical region are captured in their rest frame. We use this unique dataset, which mitigates 
complications from source variability, to intercompare the line profiles of C IV A1549, C III] A1909, Mg II A2800, 
and Ha and evaluate their implications for black hole mass estimation. We confirm that Mg II and the Balmer 
lines share similar kinematics and that they deliver mutually consistent black hole mass estimates with minimal 
internal scatter (<;0.1 dex) using the latest virial mass estimators. Although no virial mass formalism has yet been 
calibrated for C III], this line does not appear promising for such an application because of the large spread of 
its velocity width compared to lines of both higher and lower ionization; part of the discrepancy may be due to 
the difficulty of deblending C III] from its neighboring lines. The situation for C IV is complex and, because of 
the limited statistics of our small sample, inconclusive. On the one hand, slightly more than half of our sample 
(4/7) have C IV line widths that correlate reasonably well with Ha line widths, and their respective black hole 
mass estimates agree to within ~0.15 dex. The rest, on the other hand, exhibit exceptionally broad C IV profiles 
that overestimate virial masses by factors of 2-5 compared to Ha. As C IV is widely used to study black hole 
demographics at high redshifts, we urgently need to revisit our analysis with a larger sample. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — quasars: emission lines — quasars: 
general 



1. MOTIVATION 

Ever since their discovery in the early 1960s, quasars have at- 
tracted attention not only as laboratories for exploring extreme 
regimes in astrophysics, but also because they serve as useful 
beacons for probing the interstellar and intergalactic medium. 
With the growing appreciation that central black holes (BHs) 
are ubiquitous and inextricably tied to galaxy formation and 
evolution (Cattaneo et al. 2009, and references therein), quasars 
have gained even greater prominence in their unique role as 
markers of vigorous BH growth out to the highest accessible 
redshifts (Mortlock et al. 201 1; Treister et al. 201 1). 

A key development comes from our ability to estimate BH 
masses for broad-lined (type 1) active galactic nuclei (AGNs), 
including quasars, using simple parameters that can be ex- 
tracted from their rest-frame ultraviolet (UV) and optical spec- 
tra. The basic premise is that the BH mass can be approximated 
by the virial product Mbh = /^blrAV 2 /G, where /?blr, the 
radius of the broad-line region (BLR), can be estimated from 
the radius-luminosity relation calibrated through reverberating 
mapping experiments (Kaspi et al. 2000, 2005; Bentz et al. 
2009), AV is the velocity dispersion of the line-emitting gas 
that can be measured from the widths of the broad emission 
lines, G is the gravitational constant, and / is a geometric factor 
of order unity that accounts for the poorly constrained geome- 
try and kinematics of the BLR. A variety of approaches have 



been employed to calibrate the virial method for deriving Mbh 
using single-epoch spectra, using different emission lines for 
different redshift regimes. H/3 is normally the line of choice for 
low-redshift objects (Kaspi et al. 2000), as it is the principal 
line for which most reverberation mapping experiments have 
been done to date, although under some circumstances Ha is 
preferable to H/3 (Greene & Ho 2005b). At intermediate red- 
shifts, 0.75 < z < 2, McLure & Jarvis (2002) introduced a for- 
malism based on Mg II A2800, while C IV A 1549 is the only 
viable option for quasars at z > 2 (Vestergaard 2002). With 
the increasing availability of near-infrared (NIR) spectroscopy, 
these restrictions can now be circumvented by observing the 
rest-frame optical lines out to high redshift, thereby minimizing 
the additional uncertainties incurred through the extra layers of 
intermediary cross-calibrations. Nevertheless, access to NIR 
spectroscopy is still far from routine, and all extant, large spec- 
tral databases used for statistical analyses of AGNs continue to 
rely on optical surveys. 

A number of studies have investigated the robustness of dif- 
ferent broad emission lines commonly used to estimate BH 
masses. The general consensus is that Mg II serves as a reason- 
ably effective substitute for H/3, the local "standard." H/3 and 
Mg II, both low-ionization lines, are thought to arise from gas 
with common physical conditions and presumably from similar 
locations with similar kinematics. This expectation appears 
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TABLE 1 



Sample and Observations 



Name 


E{B - V) 


z 


D L 


r 


Exposure Time 


Airmass 


Seeing 




(mag) 




(Gpc) 


(mag) 


(s) 




(") 


SDSS J085615.62-001352.1 


0.035 


1.3906 


10.095 


18.88 


1800 


1.32 


0.83 


SDSS J102900.49+081325.8 


0.029 


1.3715 


9.922 


18.73 


1800 


1.41 


1.04 


SDSS J105023.68-010555.5 


0.049 


1.5349 


11.423 


19.15 


1800 


1.18 


0.77 


SDSS J122925.95-030702.4 


0.033 


1.4428 


10.572 


17.50 


3600 


1.21 


1.01 


SDSS J123816.84-005656.5 


0.023 


1.3998 


10.179 


18.79 


1800 


1.11 


0.61 


SDSS J124255.20+024957.0 


0.028 


1.4583 


10.714 


19.19 


1800 


1.13 


0.61 


SDSS J132339.23+015452.1 


0.025 


1.4921 


11.026 


19.13 


1800 


1.12 


0.64 



NOTE.— Col. (1): Official SDSS name. Col. (2): Galactic reddening, derived from Schlegel et al. (1998). 
Col. (3): Redshift, derived from the centroid of [O III] A5007, as measured from our X-shooter spectra. Col. 
(4): Luminosity distance. Col. (5): SDSS r-band PSF magnitude. Col. (6): Total integration time. Col. (7): 
Airmass during the midpoint of the observation. Col. (8): Seeing. 



to hold to a first approximation (McLure & Dunlop 2004; Sal- 
viander et al. 2007; McGill et al. 2008; Shen et al. 2008), 
and whatever small residual differences there might be appear 
correctable with empirical prescriptions (Onken & Kollmeier 
2008) or more refined spectral analysis (Rafiee & Hall 2011). 
C IV, on the other hand, turns out to be more problematic. 
While some contend that C IV can be calibrated to deliver useful 
mass estimates (Vestergaard 2002; Warner et al. 2003; Vester- 
gaard & Peterson 2006; Kelly & Bechtold 2007; Dietrich et al. 
2009; Assef et al. 201 1), others sound a more pessimistic note. 
Baskin & Laor (2005) systematically compared the profiles of 
C IV and H/3 for low-redshift (z < 0.5) quasars for which they 
could locate published optical and archival space-based UV 
spectra. The two lines agree poorly. Not only does the width 
of C IV show large and apparently systematic deviations from 
H/3, but, as long known (Gaskell 1982; Tytler & Fan 1991), 
the profile of C IV is often highly blueshifted and asymmet- 
ric, casting serious doubt as to whether the line properly traces 
gravitationally bound gafl Netzer et al. (2007), Sulentic et al. 
(2007), and Shen & Liu (2012) arrive at a similar conclusion 
from analysis of moderate-redshift (z ~ 2) quasars for which 
they secured rest-frame Balmer line measurements using near- 
IR spectroscopy; the widths of the Balmer lines exhibit little, if 
any, correlation with the widths of C IV. 

AGNs vary. Larger variations usually occur at shorter 
wavelengths and in lines of higher ionization. An impor- 
tant limitation of most of the previous studies comes from 
the fact that the rest-frame UV and optical lines were ob- 
served non-simultaneously, often separated widely apart in time 
(timescales of months to years), making any comparison be- 
tween them inherently uncertain. The width of C IV in quasars, 
for instance, changes up to ^30% on timescales of weeks to 
months (Wilhite et al. 2006a). Fortunately, variability has only 
a relatively minor impact on BH masses derived from single- 
epoch spectra. According to Wilhite et al. (2006b), Denney 
et al. (2009), and Park et al. (2012) variability affects the mass 
estimates only at the level of ^0. 1 dex, which is small, but sig- 



nificant, compared to the 0.3-0.4 dex error budget that still 
plagues the various mass estimators (Vestergaard & Peterson 
2006; McGill et al. 2008). Still, to achieve a better understand- 
ing of how the different lines most commonly used for BH mass 
estimation compare to each other, it would be highly desirable 
to perform a comparative analysis using a set of simultaneous 
rest-frame UV-optical spectra of quasars. This is the main goal 
of this paper. 




FIG. 1. — Distribution of the continuum luminosity at 3000 A and the 
FWHM of the broad component of Mg II A2800 for 1 1,015 SDSS DR7 quasars 
with z = 1.4- 1.6 (Shen et al. 2011). The objects observed with X-shooter are 
plotted as red points. The values of FWHM from SDSS were reduced by 0.05 
dex to account for the systematic difference in the method used to fit Mg II 
(Shen etal. 2011). 



8 Vestergaard & Peterson (2006) reassessed Baskin & Laor's analysis and concluded that the mismatch between C IV and H/3, though certainly substantial, is not 
as severe has had been claimed. 
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FIG. 2. — X-shooter spectra of our sample, shifted to the rest-frame of the objects. We only show the spectral regions that contain the strongest emission lines 
used for the profile analysis. The individual objects have been shifted by arbitrary constant additive offsets to improve the clarify of the plot. 



Distance-dependent quantities are calculated assuming the 
following cosmological parameters: Hq = 100 h = 71 km 
s" 1 Mpc" 1 , n m = 0.26, and n A = 0.74 (Komatsu et al. 2009). 

2. OBSERVATIONS AND DATA REDUCTION 

We selected seven relatively bright (r w 18-19 mag) quasars 
from the Seventh Data Release of the Sloan Digital Sky Sur- 
vey (SDSS DR7; Abazajian et al. 2009). As we are interested 
in simultaneous coverage of C IV A 1549 to Ha, we chose the 
targets to have z ~ 1.5, During the selection process, we in- 
spected the existing SDSS spectra to ensure that the objects 
have unambiguous broad emission lines suitable for estimating 
BH masses. We avoided sources with obvious broad absorption 
features, but aside from this, we did not apply any other selec- 
tion criteria. Figure 1 compares our sample with the distribution 
of Mg II line widths (FWHM) and 3000 A continuum luminos- 
ity for 1 1,015 z = 1 4- 1 .6 quasars contained in the DR7 catalog 
of Shen et al. (2011). Our objects are on average slightly less 
luminous than the peak of the DR7 distribution, but there are no 
obvious biases. The one apparent outlier with relatively narrow 
lines is SDSS J122925.95-030702.4, which was intentionally 
included for comparison because it has characteristics similar 
to narrow-line Seyfert 1 (NLS1) galaxies. Thus, even though 
our sample is small and somewhat ill-defined, it nonetheless 
represents an unbiased subset of luminous (L3000 > 10 45 4 erg 
s" 1 ), moderate-redshift, optically selected quasars. 

The observations were made using X-shooter (Vernet et al. 
2011), a three-arm, single-object echelle spectrograph that 



started operations in 2009 October on the VLT. The instru- 
ment covers simultaneously the wavelength range from ^3000 
A to 2.40 /im, in three arms: UVB (AA = 3000-5500 A), VIS 
(AA = 5500 A - 1.02 /mi), and NIR (AA = 1.02-2.40 /im). 
The data were taken within the framework of the French Guar- 
anteed Time under program 086.B-0320(A) (PI G. Ponti) and 
took place on 2010 February 18 UT. A summary of the obser- 
vations is given in Table 1 . 

For our observations we used slit widths of 1"3, l."2, and 
1 "2, respectively, for the three arms, resulting in resolving pow- 
ers of R = A/ A A w 4000, 6700, and 4300. The slits were 
aligned along the parallactic angle. The observations consist 
of four separate exposures of 450 s each, for a total of 1800 s, 
for all sources except SDSS J122925.95-030702.4, which was 
observed for 900 s per exposure for a total of 3600 s. The expo- 
sures were taken while nodding the object along the slit with an 
offset of 5" between exposures in a standard ABBA sequence. 
Every observation was preceded by an observation of an AO V 
telluric standard star at similar airmass. The night was not pho- 
tometricQ due to the presence of thin clouds. A high, variable 
humidity was present with precipitable water vapoiFI varying 
between 5 and 7.5 mm, rather high for Paranal but not excep- 
tional for this period of the year in Chile (so-called "Bolivian 
winter"). 

We processed the spectra using version 1 . 1 .0 of the X-shooter 
data reduction pipeline (Goldoni et al. 2006), which performed 
the following actions. As usual in processing nodding obser- 
vations, the raw frames taken at different positions were differ- 



9 archive . eso . org/asm/ambient-server ?site=paranal 

10 archive . eso . org/bin/qcl_cgi?action=qcl_plot_table&table=ambient_PWV 
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FIG. 3. — Spectral fitting for one of the objects in our sample, SDSS J1229— 0307. The original data, corrected for Galactic reddening, are plotted in black 
histograms. The fitted components include a featureless power law (dashed blue), the Balmer continuum (long dashed turquoise), and iron emission (dotted purple), 
which altogether constitute the pseudocontinuum (dot-dashed green). The final model, including fits to the principal emission lines, is shown in solid red. 



enced with each other to obtain two frames (i.e. A-B and B-A), method developed by van Dokkum (2001). The frames were 
on which cosmic ray hits were detected and corrected using the then divided by a master flat field obtained using day-time fiat 
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Parameter Units J0856-0013 J1029+0813 J1050-0105 J1229-0307 J1238-0056 J1242+0249 J1323+0154 



(9(1350 - 2200 A) 






-1.37 ±0.11 


/3(2200 - 3500 A) 






-1.49 ±0.12 


/3(4200 - 5600 A) 






-2.50 ±0.20 


/3(5600 - 7500 A) 






-0.83 ±0.066 


log AL^(1350A) 


erg s~ 


-l 


45.75±0.022 


log AL^,(3000A) 


erg s~ 


_1 


45.37±0.022 


log AL^(5100 A) 


erg s~ 


-l 


45.25±0.022 


vc IV 


km s~ 


i 


-846±149 


FWHMc iv 


km s~ 


i 


4758±714 


log L c iv 


erg s~ 




44.29±0.072 


VC III] 


km s" 


i 


-693±342 


FWHM C m] 


km s~ 


i 


5878±1176 


log L c in] 


erg s~ 


-l 


43.83±0.043 


fMg II 


km s~ 


i 


240±49 


FWHM Mg ii 


km s~ 




3974±397 


log Lm e ii 


erg s~ 


-l 


43.75±0.043 


FWHM [0 in] 


km s~ 


i 


612±31 


l0gi[ O III] 


erg s~ 


-l 


43.31±0.043 


«Ha 


km s~ 


i 


-2±90 


FWHMhc 


km s" 


i 


3951±198 


log Lhc 


erg s" 


-l 


44.14±0.021 


log Fe II A4570/H/3 






-0.12±0.055 


log Fe II UV/Mg II 






0.50±0.063 


logM BH (C IV) 


M & 




8.94±0.4 


log Mbh (Mg II) 


Me 




8.84±0.3 


logMBH(Ha) 


Me 




8.85±0.3 


lOgibol/iEdd 






-0.71±0.3 



-1.14 ±0.091 


-0.73 ± 0.058 


-1.55 ±0.12 


-1.91 ±0.15 


-2.48 ±0.19 


-1.80 ±0.14 


-1.55 ±0.12 


-1.54 ±0.12 


-2.25 ±0.18 


-1.45 ±0.12 


-1.06 ±0.085 


-1.53 ±0.12 


45.75±0.022 


45.69±0.022 


46.32±0.022 


45.66±0.022 


45.55±0.022 


46.05±0.022 


45.44±0.022 


45.55±0.022 


45.98±0.022 


154±180 


842±167 


-1102±497 


3233±485 


4193±629 


5570±836 


44.67±0.072 


43.75±0.072 


44.59±0.072 


-207±202 


147±737 


-273±887 


2625±525 


6121±1224 


5931±1186 


43.77±0.043 


43.89±0.043 


44.25±0.043 


303±202 


551±362 


482±286 


2995±300 


2763±276 


2262±226 


43.77±0.043 


43.65±0.043 


43.99±0.043 


560±28 


558±28 


1300±65 


43.53±0.043 


43.11±0.043 


43.51±0.043 


61±24 


260±147 


575±211 


3058±153 


3891±195 


2479±124 


44.41±0.021 


44.26±0.021 


44.64±0.021 


-0.80±0.055 


-0.37±0.055 


-0.12±0.055 


0.50±0.063 


0.93±0.063 


0.66±0.063 


8.61±0.4 


8.80±0.4 


9.38±0.4 


8.74±0.3 


8.62±0.3 


8.69±0.3 


8.74±0.3 


8.63±0.3 


8.66±0.3 


-0.41±0.3 


-1.13±0.3 


0.21±0.3 



-2.07 ± 0.17 


-1.67 ±0.13 


-1.63 ±0.13 


-1.77 ± 0.14 


-2.15 ±0.17 


-1.42 ±0.11 


-2.36 ±0.19 


-2.09 ±0.17 


-1.89 ±0.15 


-1.07 ±0.086 


-0.63 ±0.050 


-1.19 ±0.095 


45.96±0.022 


45.93±0.022 


45.80±0.022 


45.55±0.022 


45.53±0.022 


45.58±0.022 


45.46±0.022 


45.37±0.022 


45.35±0.022 


-795±601 


-700±230 


-689±459 


6338±951 


5875±881 


3638±546 


44.22±0.072 


44.49±0.072 


44.20±0.072 


-432±294 


324±601 


-482±134 


5244±1049 


6230±1246 


2329±466 


43.73±0.043 


43.68±0.043 


43.73±0.043 


163±207 


26±237 


-184±158 


3802±380 


4421±442 


3301±330 


43.73±0.043 


43.86±0.043 


43.61±0.043 


630±32 


1260±63 


450±23 


43.15±0.043 


42.99±0.043 


42.87±0.043 


211±182 


45±38 


-337±342 


3446±172 


3945±172 


3871±194 


44.21±0.021 


44.19±0.021 


44.13±0.021 


-0.31±0.055 


-0.18±0.055 


-0.41±0.055 


0.56±0.063 


0.45±0.063 


0.48±0.063 


9.30±0.4 


9.21±0.4 


8.74±0.4 


8.89±0.3 


9.02±0.3 


8.78±0.3 


8.76±0.3 


8.87±0.3 


8.83±0.3 


-0.41±0.3 


-0.61±0.3 


-0.59±0.3 



NOTE. — j3 — index of the power-law continuum in the wavelength ranges shown. XL\ = monochromatic luminosity at the given wavelengths, v — 
velocity shift of the line, defined by the peak of the line relative to the systemic velocity as set by [O III] A5007. FWHM — full width at half maximum 
measured from the best-fit model of the line profile. L — luminosity of the emission line. Fe II A4570/H/5 — luminosity ratio of the optical Fe II A4570 
feature (4434-4684 A) to broad H/3. Fe II UV/Mg II = luminosity ratio of the UV Fe II feature (2200-3090 A) to broad Mg II. Mbh = BH mass using 
the C IV, Mg II, and He* virial mass estimators, as explained in Section 4. Lboi/^Edd — Eddington ratio, estimated assuming Lboi — 9.8ALa(5100 A). 



field exposures with halogen lamps. The orders were extracted 
and rectified in wavelength space using a wavelength solution 
previously obtained from calibration frames. The resulting rec- 
tified orders were then shifted and added to superpose them, 
thus obtaining the final two-dimensional spectrum. The orders 
were then merged, and, in the overlapping regions, the merg- 
ing was weighted by the errors that were propagated during the 
process. From the resulting two-dimensional, merged spectrum 
a one-dimensional spectrum was extracted at the source's po- 
sition. The one-dimensional spectrum with the corresponding 
error file and bad pixel map is the final product of the reduction. 

To perform flux calibration, we used different procedures for 
the UVB data and for the VIS-NIR data. In the UVB band we 
used an observation of the flux standard GD 7 1 (Bohlin 2007) 
taken in the beginning of the night. We reduced the data us- 
ing the same steps as above, but in this case we subtracted the 
sky emission lines using the method of Kelson (2003). The ex- 
tracted spectrum was divided by the flux table of the star from 
the CALSPEC HST databasfj to produce the response func- 
tion, which was then applied to the spectrum of the science tar- 
gets. For the VIS and NIR arms, we used A0V stars as both 
flux and telluric standards. We extracted the A0V spectra with 
the same procedure used for the flux standard and used these 
spectra to apply telluric corrections and flux calibrations simul- 
taneously, using the package Spextool (Vacca et al. 2003). 
We verified that the final spectra of the three arms were com- 
patible in the common wavelength regions and then adjusted 

11 www .stsci.edu/hst/ observatory/ cdbs/calspec. html 



the mean continuum flux in the overlap range between 4000 
A and 9000 A to be consistent with the SDSS spectra, which 
have a more reliable absolute flux calibration. Figure 2 shows 
the calibrated, rest-frame spectra of the quasars observed in our 
program. 

3. SPECTRAL FITTING 

We correct the X-shooter and SDSS spectra for Galactic 
extinction using the extinction map of Schlegel et al. (1998) 
and the reddening curve of Fitzpatrick (1999). The spectra 
are transformed into the rest frame using the redshift as de- 
termined from the peak of the best-fit model (see below) for 
the [O III] A5007 line. To improve on the absolute flux calibra- 
tion, we rescale the X-shooter spectra to the flux density level 
of the SDSS spectra over the wavelength region where the two 
overlap. Here we present a brief description of the spectral fit- 
ting, which is based on ^-minimization using the Levenberg- 
Marquardt technique within the MPFIT package (Markwardt 
2009). 

The analysis for the optical spectra closely follows the 
methodology for decomposition of AGN spectra described in 
Dong et al. (2008). We do not correct for starlight contam- 
ination, which is expected to be negligible for quasars in our 
luminosity range (L 30 oo > 10 45 - 4 erg s" 1 ). As the broad emis- 
sion lines, particularly the Fe II multiplets, are so broad and 
strong that they merge together and essentially leave no line- 
free wavelength regions, we fit simultaneously the nuclear 
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FIG. 4.— Model line profiles, shown on a velocity scale, for Ha A6563 (black, solid), Mg II A2800 (red, dotted), C III] A1909 (green, dashed), and C IV A1549 
(blue, long dashed). The zero point of the velocity scale is referenced with respect to [O III] A5007. 



continuum, the Fe II multiplets, and other emission lines. The 
featureless continuum of type 1 AGNs is not well described by 
a single power law when a large range of wavelengths is con- 
sidered (e.g., Vanden Berk et al. 2001). We approximate it by a 
broken power law, with free indices for the Ha and H(3 regions. 
The optical Fe II emission is modeled with two separate sets 
of analytical spectral templates, one for the broad-line system 
and the other for the narrow-line system, constructed from mea- 
surements of IZw 1 by Veron-Cetty et al. (2004). Within each 
system, the respective set of Fe II lines is assumed to have no 
relative velocity shifts and the same relative strengths as those 
in I Zw 1. Emission lines are modeled as multiple Gaussians. 
Following Dong et al. (2011), we assume that the broad Fe II 
lines have the same profile as broad Ha. In our data set, the 
Hf3 region is significantly noisier than the Ha region. In our 
subsequent analysis we will use Ha in lieu of H/3; Greene & 
Ho (2005b) have shown that the broad component of Ha is on 



average slightly narrower than H/3, but the two closely track 
each other. All narrow emission lines are fitted with a single 
Gaussian, except for the [O III] AA4959,5007 doublet, each of 
which is modeled with two Gaussians, one accounting for the 
line core and the other for a possible blue wing seen in many 
objects (Greene & Ho 2005a). 

For the near-UV spectra, we focus our analysis on the region 
around Mg II AA2796, 2803, following the method described in 
Wang et al. (2009). Here the pseudocontinuum consists not 
only of a local power-law continuum and an Fe II template, but 
also an additional component for the Balmer continuum. The 
Fe II emission is modeled with the semi-empirical template for 
I Zw 1 generated by Tsuzuki et al. (2006). To match the width 
and possible velocity shift of the Fe II lines, we convolve the 
template with a Gaussian and shift it in velocity space. As in 
Dietrich et al. (2002), the Balmer continuum is assumed to be 
produced in partially optically thick clouds with a uniform 
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FIG. 5. — Comparison of FWHM of Ha with FWHM of (a) Mg II, (b) C IV, and (c) C III]. The error bars represent measurement uncertainties as given in Table 
2 (see Section 3). 



temperature. Each of the lines of the Mg II doublet is modeled 
with two components, one broad and the other narrow. The 
broad component is fit with a truncated, five -parameter Gauss- 
Hermite series; a single Gaussian is used for the narrow com- 
ponent. The two doublet lines are assumed to have identical 
profiles, a fixed separation set to the laboratory value, and a 
flux ratio A2796/A2803 set to be between 2:1 and 1:1 (Laor et 
al. 1997). 

The fitting for the UV region is performed with a modified 
version of a code initially written and kindly provided by Jian- 
Guo Wang. We employ the UV Fe II + Fe III template for 
I Zw 1 generated by Vestergaard & Wilkes (2001). In light 
of the moderate signal-to-noise ratio of our data, we adopt the 
same scaling factor for the Fe II and Fe III emission; this as- 
sumption seems adequate from visual inspection of the fits. As 
in Vestergaard & Wilkes (2001), the power-law continuum and 
iron emission (both constituting the pseudocontinuum) are fit to 
the emission-line-free spectral regions from 1300 A to 2200 A. 
After the pseudocontinuum is subtracted, we concentrate our 
fits on two regions: (1) — 1350- 1700 A, which contains the 
blend of Si IV and O IV] at -1400 A, C IV A 1549, He II A 1640, 
OIII] A 1663, and Al II A 1670, and (2) - 1800-2000 A, which 
covers Si IIA1817, Al III AA1855, 1863, Si III] A1892, C III] 
A1909, and the Felll UV34 triplet at -1914 A. We model the 
emission lines with multiple Gaussian components, using the 
minimum number necessary to achieve a satisfactory fit within 
the signal-to-noise constraints of the data. The C IV profile is 
fit with three Gaussians, taking care to avoid narrow absorption 
features when present (as in J1238-0056 and J1323+0154; see 
Figure 2). For the C III] region, Si II, Al III, Si III], and Felll 
are each fit with a single Gaussian, while C III] itself is fit with 
two Gaussians, which in most cases yield better residuals than 
a single Gaussian. This suffices for our purposes. Our primary 
objective is to obtain a robust characterization of the profile of 
C IV and C III], not to achieve a detailed model for every line in 
this tremendously complicated spectral region. 

Figure 3 illustrates the fits for one of objects, and the results 
of the fits for the entire sample are summarized in Table 2. The 
emission-line luminosities, velocity shifts (peak of the line; v), 
and velocity widths (FWHMs) are measured from the best-fit 
models of the line profiles, except for the case of Mg II, whose 
line peak and FWHM are measured from the model of the sin- 
gle doublet line, Mg II A2796. The strength of the optical 



Fe II A4570 emission is integrated over the wavelength range 
4434-4684 A, and that for U V Fe II is integrated over the range 
2200-3090 A. For all measured emission-line fluxes, we regard 
the values as reliable detections when they have greater than 3 a 
significance. As the instrumental resolution (FWHM m 45-75 
km s" 1 ) is small compared to the measured line widths (FWHM 
s=s 500-6000 km s -1 ), no correction for instrumental resolution 
has been applied to the line widths. 

We follow the bootstrap method of Dong et al. (2008; their 
Section 2.5) to estimate measurement uncertainties. The typ- 
ical 1 a errors on the fluxes of the strong lines considered in 
this paper are quite small, <;10% (see Dong et al. 2008, 2011; 
Wang et al. 2009). We adopt 10% errors for the fluxes of Mg II, 
C III], and [O III] and 5% for Ha. The only exception is C IV, 
whose final flux depends on the exact procedure adopted to fit 
the "red shelf" (region — 1600- 1700 A; see, e.g., Fine et al. 
2010). We assume that the red shelf is intrinsic to C IV; ex- 
cluding it reduces the inferred line flux by —10-25%. To be 
conservative, we adopt an error of 18% for the flux of C IV. 
The FWHM values are best measured for Ha (error 5%), fol- 
lowed by Mg II (error 10%), C IV (error 15%), and C III] (error 
15%-25%). C III] is particularly problematic because of se- 
vere blending with Al III A A 1 855 , 1 863, Si III] A 1 892, and Fe III 
UV34. 

4. RESULTS 

Figure 4 shows the final, fitted profiles for the sample, nor- 
malized by their amplitude and plotted on a velocity scale. A 
direct comparison among the line widths is given in Figure 5. 
Several trends are immediately obvious. 

1. Ha shows, by far, the most symmetric profiles and the 
least amount of systemic velocity shift. This is not un- 
expected, given the overall similarity between Ha and 
H/3 profiles (Greene & Ho 2005b), reinforcing the com- 
monly adopted view that the Balmer lines offer the most 
reliable BH virial mass estimators (but see Vestergaard 
et al. 201 1 for a different point of view). 

2. Consistent with other studies (e.g., McLure & Dunlop 
2004; Shen et al. 2008; Shen & Liu 2012), we find that 
Mg II on average serves as a reasonably decent proxy for 
the Balmer lines, at least for five out of the seven objects 
in our sample. The two exceptions are J1050-0105 and 
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J1229-0307, whose Mg II lines are distinctly narrower 
than Ha and somewhat asymmetric. 

3. Among the five objects that show better match be- 
tween Mg II and Ha, the lines agree best near the core 
(FWHM) but often diverge noticeably in the wings. Ha 
generally has more extended wings than Mg II. This im- 
plies that comparison between virial mass estimators for 
Mg II and the Balmer lines depends sensitively on the 
adopted measure of line width. The agreement should 
be reasonably close for the FWHM (Figure 5(a); mean 
and standard deviation (FWHM(Mg II)/FWHM(Ha)) = 
0.95 ±0.14 km s" 1 ), but the line dispersion, more sensi- 
tive to the line wings, will be systematically higher for 
the Balmer lines compared to Mg II. 

4. The behavior of C IV is much more perplexing. Apart 
from the presence of a prominent extended red wing — 
a model-dependent feature sensitive to the adopted fit- 
ting procedure (Fine et al. 2010) — the line peak ex- 
hibits a range of velocity shifts, from ^800 km s" 1 
to -1100 km s . A direct comparison between C IV 
and Ha line widths reveals two clusters of points (Fig- 
ure 5(b)). Four of the quasars have C IV line widths 
that are roughly consistent with those of Ha (mean 
and standard deviation (FWHM(C IV)/FWHM(Ha)) = 
1.07 ± 0.11 km s" 1 ), whereas the other three 
(J1229-0307, J1238-0056, and J 1242+0249) are signif- 
icantly broader, (FWHM(C IV)/FWHM(Ha)) = 1.86 ± 
0.38 km s" 1 . Both are difficult to understand in 
the context of a radially stratified BLR with a ve- 
locity field dominated by Keplerian rotation, wherein 

1 /2 

AV oc R^ R . Not many AGNs have been reverbera- 
tion mapped in C IV, but the handful for which ade- 
quate data exist show that C IV arises from a region 
that is a factor of ~ 2 more compact than H/3 (e.g., 
Korista et al. 1995; Onken & Peterson 2002). If 
these results apply to higher luminosity quasars, they 
imply that FWHM(C IV) « \/2FWHM(Ha) [for sim- 
plicity, we assume FWHM(Ha) = FWHM(H/3)]. Our 
sample, albeit small, does not seem to conform to this 
simple expectation. The three outliers with unusu- 

12 We assume bolometric luminosity Lbol = 9.8AL^(5100) (McLure & Dunlop 
Mbh estimated from Ha. 



ally high FWHM(C IV)/FWHM(Ha) share one com- 
mon characteristic: the peak of the C IV is signifi- 
cantly blueshifted, by -700 to -1100 km s" 1 . How- 
ever, not all objects with C IV blueshifts show enhanced 
FWHM(C IV)/FWHM(Ha); two of the other four ob- 
jects also show blueshifts of comparable magnitude. 
Nor is excess C IV width uniquely associated with any 
obvious AGN property. While J 1229-03 07, with its 
relatively narrow lines [FWHM(Ha) « 2500 km s" 1 ], 
strong Fe II emission (log Fe II A4570/H/3 = -0.12), and 
high Eddington raticpl (log Lboi/^Edd = 0.21) qualify it as 
a NLS1, whose C IV line may be particularly problem- 
atic for mass determination (Vestergaard et al. 2011), 
the other two outliers do not. 

5. C III] poses an even greater challenge to understand. 
While there are no large velocity shifts, the width of 
C III] appears to be completely erratic (Figure 5(c)). 
Two of the objects have unusually narrow cores that 
result in FWHM(C III])/FWHM(Ha) < 1, and the rest 
are characterized by FWHM(C III])/FWHM(Ha) > 1. 
Counterintuitively, FWHM(C III]) > FWHM(C IV) in 
four out of the seven sources. This agrees with the re- 
sults of Brotherton et al. (1994), Jiang et al. (2007), 
and Greene et al. (2010), but is inconsistent with 
the study of Shang et al. (2007), who generally find 
FWHM(C III]) ^ FWHM(C IV). Using a much larger 
sample, Shen & Liu (2012) find no obvious offset be- 
tween FWHM(C III]) and FWHM(C IV), although the 
two are poorly correlated. The discrepancy associated 
with C III] may arise, at least in part, from the uncer- 
tainty in deblending the line from its surrounding con- 
taminating features. 

How do these results impact BH mass determinations? Fig- 
ure 6 graphically illustrates the answer. We estimate BH masses 
using the C IV-based formalism of Vestergaard & Peterson 
(2006), the Mg Il-based formalism of Vestergaard & Osmer 
(2009), and the Ha-based formalism of Greene & Ho (2005b), 
as updated by Xiao et al. (201 1) to account for the latest BLR 
size-luminosity relation of Bentz et al. (2009). For consistency 
with our adopted C IV and Mg II mass estimators, we further 

■), Eddington luminosity Z-Edd = 1.26 X 10 38 (Mbh/Mq) erg s _1 , and BH mass 
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adjust Xiao et al.'s prescription so that its virial coefficient / 
matches the value advocated by Onken et al. (2004). We do 
not consider C III] because there is currently no mass estimator 
based on this line, and we feel discouraged by the line width 
comparison presented above. The BH masses derived from 
these virial mass estimators typical have statistical uncertain- 
ties of ~ 0.3 - 0.4 dex. For concreteness, we adopt an error bar 
of 0.4 dex for M BH (C IV) and 0.3 dex for both M BH (Mg II) and 
M BH (Ha). 

Vestergaard et al. (201 1) stress the importance of using mass 
estimators calibrated on the same mass scale for a proper com- 
parison between masses derived from different lines. Both the 
C IV and Mg II masses are tied to a common scale based on 
reverberation-mapped AGNs. Vestergaard & Osmer, however, 
model the broad Mg II line as a single component, whereas we 
follow Wang et al. (2009) and treat Mg II as a doublet, which re- 
sults in slightly narrower line widths. These two methods yield 
Mg II line widths that differ on average only by <~ 0.05 dex 
(Shen et al. 201 1); we scale up all of our Mg II line widths by 
this constant factor to account for the minor systematic offset. 
To date the Ha BH mass estimator of Greene & Ho (2005b) has 
not been calibrated directly against the reverberation-mapped 
AGNs. Fortunately, the analysis of Shen et al. (2011) indi- 
cates that the Ha-based masses show only a mean offset of 
0.08 dex with respect to Vestergaard & Peterson's H/3-based 
masses, which, like the C IV and Mg II masses, are calibrated to 
the same scale tied to the reverberation-mapped AGNs. For the 
purposes of this paper, we will not worry about this small dis- 
crepancy, which does not impact any of our main conclusions. 

As foreshadowed by the line width analysis, Ha and Mg II 
deliver reasonably consistent mass estimates, which for our 
sample spans a small range around M B h ~ 10 87±a2 Mq. For 
our choice of mass estimators, M B n(Mg II) agrees very well 
with M BH (Ha): (log M BH (Mg II)-log M BH (Ha)) w 0.03 ± 0.07 
dex. By contrast, the comparison between M BH (C IV) and 
M B h(Ho:) is less clear-cut. The three objects (labeled in Fig- 
ure 6) with anomalously broad, blueshifted C IV profiles all 
have M BH (C IV) in excess of 10 9 Mq, deviating fromM BH (Ha) 
by 0.34 dex to as much as 0.72 dex. The most discrepant ob- 
ject is the NLS1 J1229-0307, confirming Vestergaard et al.'s 
(201 1) suspicion that C IV masses are especially unreliable for 
this class of AGNs, but the other two show no obvious warn- 
ing signs as to why C IV should misbehave. All three C IV 
outliers appear relatively normal in the Mg II versus Ha com- 
parison (Figure 6(a)). The remaining four objects fare better: 
(log M BH (C IV)-log M BH (Ha)) w 0.01 ± 0.14 dex. While this 
may be regarded as reassuring confirmation that C IV-based 
masses can be trusted in at least some objects, the difficulty is 
that, in the absence of independent evidence from lower ioniza- 
tion lines (Mg II, H/3, or Ha), we have no means of forecasting 
which objects are reliable or not. This result is unsettling, espe- 
cially in light of the very large masses (M BH ~ 10 9 - 10 10 Mq) 
routinely inferred for high-redshift quasars and the astrophysi- 
cal implications attached to them. 

5. SUMMARY 

The widths of broad emission lines in active galaxies, when 
combined with physical dimensions inferred from the size- 
luminosity relation empirically calibrated from reverberation 
mapping experiments, provide a powerful and efficient means 
of estimating BH masses for large samples of sources detected 
at all redshifts. Over the past decade, a number of virial BH 
mass estimators have been devised and extensively used to de- 



termine BH masses, facilitating a wide range of investigations 
on BH demographics and AGN physics. Despite their popu- 
larity, however, there are nagging doubts as to the reliability 
of these mass estimators. While calibrations based on low- 
ionization lines, especially Ha and H/3, are reasonably secure, 
higher redshift observations often depend on UV lines that are 
less well understood. The most contentious mass estimator is 
that based on C IV A 1549, the workhorse for studies of high- 
redshift quasars, as the kinematics of this line may be signif- 
icantly affected by winds and other non- virial motions. Mg II 
appears to be safer, but even it may not be completely immune 
to systematic biases. 

Spectral variability complicates the comparison between dif- 
ferent mass estimators, if they derive from spectra taken at 
different times. To mitigate this effect, we have undertaken 
an experiment to acquire spectra of a small sample of seven 
moderate-redshift quasars that simultaneously cover the rest- 
frame U V through optical spectral regions ( 1 300-7500 A). This 
dataset is enabled by the unique capabilities of the X-shooter in- 
strument on the VLT, which delivers simultaneous spectra from 
^3000 A to 2.5 /im. At z ~ 1.5, this allows us to access the 
principal broad emission lines from C IV to Ha. 

In accord with other studies, we find that Mg II and the 
Balmer lines (this study uses Ha instead of H/3) have similar 
velocity widths near the core (FWHM) of their profiles, but 
Ha generally has more extended, higher velocity wings than 
Mg II. Mg Il-based and Ha-based BH masses agree to better 
than ~ 0.1 dex. The C III] A1909 line widths are difficult to 
interpret. Contrary to naive expectations, FWHM(C III]) > 
FWHM(C IV) in most of our objects, but some also have un- 
usually narrow lines, narrower than even those of Ha. While 
these discrepancies can perhaps be attributed to the difficulties 
of line deblending in the highly crowded spectral region near 
^1900 A, it appears that it would be challenging to devise a 
robust BH virial mass estimator using C III]. The verdict on 
C IV is mixed. While roughly half of the C IV-based masses are 
reasonably consistent with those derived from Ha, the others 
are systematically high by factors of 2-5. These extreme out- 
liers not only have unusually broad C IV line widths, but their 
line peaks are all systematically blueshifted by several hundred 
to a thousand km s" 1 , suggesting that a significant fraction of 
the emission arises from outflowing or dynamically unrelaxed 
gas. However, systemic C IV blueshifts are a common feature 
in AGNs, and they do not appear to be a clean predictor of 
which objects show deviant C IV line widths or masses; two of 
the objects whose C IV-based masses agree with those derived 
from Ha also have systematic C IV blueshifts of comparable 
magnitude. We are thus left in an uncomfortable predicament: 
in the absence of independent confirmation from lower ioniza- 
tion lines, we do not know, a priori, which C IV profiles provide 
more accurate mass estimates. 

To end on a more positive note, we emphasize that our results 
clearly suffer from small-number statistics. Other investigators 
(e.g., Vestergaard & Peterson 2006; Greene et al. 2010; As- 
sef et al. 2011) are more optimistic that C IV can be used to 
study BH demographics at high redshifts. It would be impor- 
tant to secure simultaneous, or at least near-contemporaneous, 
rest-frame U V-optical spectra for a larger, better defined s am- 
ple of AGNs to revisit the issues raised in this study. Be- 
cause high-luminosity AGNs vary on time scales of weeks to 
months, complete simultaneity in spectral coverage (such as 
those presented here) is desirable, but not truly necessary. Near- 
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contemporaneous observations, taken within a span of a few 
days (e.g., during the same observing run or coordinated be- 
tween two different telescopes), would suffice. 
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